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I.  INTRODUCTION 


Thermal  Insulation  characterization  of  polymers  is  an  Important  require¬ 
ment  for  present  and  future  aerospace  systems.  Examples  Include  a  Neoprene/ 
Kevlar  composite  used  in  rocket  motors,  as  well  as  insulating  blankets  of 
fused  silica  or  other  materials.  Furthermore,  the  construction  of  new  homes 
requires  the  use  of  lightweight  fibrous  insulation  whose  thermal  properties 
need  to  be  known. 

The  significance  of  thermal  radiation  in  fibrous  Insulations  has  been 
I  2 

demonstrated  •  by  conducting  guarded  hot-plate  experiments  under  vacuum. 

Tong  and  Tien^  categorized  the  analyses  of  such  experiments  as  those  which 
experimentally  determine  effective  conductivities  caused  by  radiation  and 
those  which  consider  the  detailed  equations  governing  the  intensity  of 
radiation  in  an  absorbing  and  scattering  medium. 

The  thermal  diffusivity  a  is  the  parameter  of  interest  in  transient-heat 
conduction  problems,  as  opposed  to  steady-state  problems  in  which  only  the 
thermal  conductivity  K  plays  a  role.  For  a  homogeneous  body,  a  may  be  calcu¬ 
lated  from  K,  the  density  p,  and  the  specific  heat  c.  If  a,  p,  and  c  are 
measured  in  a  given  temperature  range,  K  may  be  calculated  over  that  range. 
This  is  especially  advantageous  when  steady-state  calorimetry  for  K  is 
difficult,  e.g.  at  high  temperatures. 

One  method  of  measuring  the  thermal  diffusivity  is  the  flash  tech- 

S—  1 2 

nique  .  A  flash  of  radiant  energy  is  deposited  over  one  face  of  a  disk¬ 
shaped  sample  that  is  surrounded  by  a  heating  element  and  enclosed  in  a 
vacuum.  For  a  homogeneous  sample,  a  may  be  calculated  from  the  sample  thick¬ 
ness  L  and  the  time  tii2  (the  time  for  the  rear  face  to  achieve  one  half  its 
ultimate  temperature).  To  within  a  numerical  factor,  a  is  essentially 
L^/tj/z. 

12 

It  has  been  observed  that  problems  involving  large  temperature 
differences  and  short  heat-pulse  propagation  times  may  lead  to  situations  in 
which  the  basic  equations,  boundary  conditions,  and  assumptions  of  the  flash 


technique  data-reductlon  scheme  are  inapplicable.  For  these  problems  there 
may  be  no  way  to  convert  an  observed  oscilloscope  trace  into  a  meaningful 
value  for  some  physical  parameter.  Conditions  occurring  during  space  vehicle 
reentry,  laser  irradiation  problems,  measurements  at  cryogenic  temperatures, 
and  materials  with  anisotropic  conductivity  may  lead  to  situations  in  which 
the  flash  dlffuslvlty  method  may  not  be  viable. 

On  the  other  hand,  some  nonhomogeneous  materials  have  been  successfully 
characterized  by  the  flash  technique.  The  flash  technique  has  been  applied  to 
some  heterogeneously  dispersed  composites  and  layered  samples.  Including  some 
carbon  fibers,  fiber-reinforced  materials.  Individual  layers  of  layered 
composites,  and  dispersed  composites. 

The  continuing  challenge  of  composite-sample  characterization  by  means 
of  the  flash  technique  Is  Illustrated  by  recent  case  histories  and  journal 
publications.  Flash-technique  analysls^^  of  a  Neoprene-covered  Kevlar  cloth 
composite  revealed  that  the  dlffuslvlty  value  could  not  be  calculated  for  the 
Kevlar  layer  solely  on  the  basis  of  the  measured  rise  times  and  the  properties 
of  the  virgin  Neoprene.  The  rise  times  were  too  short.  Indicating  that  the 
conductivity  values  for  the  virgin  Neoprene  were  too  small  for  the  Neoprene  in 
the  composite.  In  a  second  example,  fibrous  Insulation,  it  has  been  shown 
that  steady-state  and  transient  methods  of  determining  thermal  conductivity 
often  give  different  results,  especially  for  such  porous  Insulating  materials 
as  rockwool  and  kapok. ^  This  discrepancy  arises  from  the  contribution  of 
Internal  radiation  to  the  conduction  of  heat  and  should  be  even  larger  for 
an  intrinsically  transparent  material.  Although  recent  publications^*^  have 
addressed  this  Issue  from  other  viewpoints.  Internal  radiation  has  not,  to  the 
present  author's  knowledge,  previously  been  allowed  for  in  the  development  of 
the  flash  technique. 

Therefore,  we  see  that  some  currently  interesting  sample  types  are  In 
some  way  too  complicated  for  available  flash-technique  analysis.  The 
instrumentation  Is  producing  experimental  outputs  that  are  apparently  In  need 
of  more  refined  Interpretation.  In  any  particular  case,  any  one  of  a  number 
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of  different  adjustments  may  be  necessary.  These  adjustments  are  necessitated 
by  four  categories  of  physical  effects: 

1.  Radiant-energy  penetration  during  flash  heating  of  the  dlffuslvlty 
sample  (time  =0). 

2.  The  contribution  of  Internal  radiation  and  absorption  to  heat 
transfer  (time  >0). 

3.  Heterogeneous  sample  media. 

4.  Anisotropic  sample  media. 

An  adjustment  accounting  for  the  first  effect  and  applying  to  both  the 
layered-composite  and  porous-insulation  types  Is  discussed  In  this  paper. 

This  mathematical  modification  of  the  current  flash-technique  data-reduction 
method  explicitly  allows  for  an  exponentially  decaying  spatial  penetration  of 
Incident  radiant  energy  into  a  layered  sample.  A  method  of  experimentally 
evaluating  the  characteristic  length  associated  with  the  exponential 
penetration  function  will  be  proposed. 


II.  THEORETICAL  DEFINITION  OF  THE  PROBLEM 


A.  HOMOGENEOUS  SLAB 

For  a  homogeneous  slab  between  x  =  0  and  x  »  L,  the  governing  equations 
(as  in  Ref.  14)  are 


J  =  -KVT  (1) 

7  •  J  +  pc  *  0  (2) 

where  T  =  T(x,t)  is  the  temperature  as  a  function  of  position  x  and  time  t, 

K,  p,  and  c  are  the  thermal  conductivity,  density,  and  specific  heat,  respec¬ 
tively,  and  J  is  the  flux  (energy  per  area  and  time).  Homogeneity  implies 
constant  K,  p,  and  c,  so  combining  the  equations  gives  the  partial  differ¬ 
ential  equation  (PDE) 


2 

aV^T 


3t 


(3) 


where 


a 


pc 


(4) 


The  boundary  conditions  (BC)  for  the  homogeneous  slab  are  7T  =  0  at 
X  =  0  and  L  (no  heat  loss).  The  initial  conditions  (IC)  are 


T(x,0)  =  6(x)  (5) 

pc 

where  Q  is  the  energy  per  area  deposited  on  the  slab  surface  (x  *  0  at 
t  •  0),  and  6(x)  is  the  Dirac  delta  function.  Separation  of  variables 
[T(x,t)  E  W(x)T(t)]  leads  to  the  ordinary  differential  equations  (ODE) 

W"  +  (X/a)W  ■  0  and  t'  Xt  ■  0,  where  X  is  the  separation  constant  and  primes 
denote  differentiation.  In  general  for  this  Sturm-Llouville  problem. 


Wjj  “  Ajj  sin  BjjX  +  Bj^  cos  B^x  and  t  =  e  =  e  n  .  The  above  BC  imply  that 
all  “  0 ,  and  the  eigenvalue  equation  sin  B^L  =  0  leads  to 


B^  n  =  0,1.2... 


(including  negative  n  only  changes  the  B^^  definition).  The  B^^  are  specified 
by  the  IC  with  the  result 


2  2 

nirx  /-  n  ir 


[l  +  2  I  cos  exp  ( 


The  ultimate  x  =  L  temperature  is  Q/pcL,  so  the  dimensionless  temperature 
versus  time  profile  at  x  =  L  is 

V  =  1  +  2  I  exp  [— 


Approximating  the  V  result  by  the  first  two  terms  (valid  as  t  +  •) ,  one 
observes  that  V  =  1/2  when 

2  2 
L  An4  „  _  L 

a  =  — s -  =  0.1404  -  (9 


Using  the  full  series  expression. 


0.1388 


Thus,  the  thermal  diffusivity  of  a  homogeneous  slab  can  be  measured  by 

measuring  the  thickness  L  and  the  rise  time  t^/2*  This  can  be  done  to  within 

0.5%  precision. On  the  other  hand,  the  IC  and  BC  assumptions,  such  as  delta 

function  temperature-pulse  across  the  entire  sample  and  complete  thermal 

isolation,  limit  the  accuracy  of  the  measurement.  Models  have  been  developed 
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that  generate  corrections  up  to  21%  in  magnitude  in  the  measured  values. 


B.  COMPOSITE  SLAB 

Surprisingly  enough,  although  a  unique  solution  can  also  be  shown  to 
exist  for  a  composite  slab  consisting  of  individual  layers,  that  solution 
cannot  be  constructed  by  an  extension  of  the  Sturm-Llouville  approach 
(p.  A14  of  Ref.  15).  Instead,  the  Laplace  transform  technique  (Chapter  Xll 
in  Ref.  14)  is  used.  The  standard  analysis  for  the  simplest  composite  (two 
layers  and  one-dimensional  heat  flow)  will  first  be  summarized,  using  mostly 
the  notation  of  Ref.  10.  Then,  in  the  following  sections,  recommendations  for 
modifying  the  mathematics  and  the  measurement  algorithm  will  be  stated, 
relying  on  the  standard  analysis  as  a  starting  point. 

The  two-layer  problem  is  defined  by  the  following  equations: 

Partial  differential  equations  (PDE): 

9T 

KiV^Ti  -  p^Cj  <  X  <  0  (11) 

3T 

■  V2  ar  “<■'<‘2  “2) 

Initial  conditions  ( IC) : 


T^(x,0)  -  0  -jl^<x<0 

T2(x,0)  -  0  0  <  X  < 

Boundary  conditions  (BC); 

♦ 

-K^7Tj(-£^ ,t)  ■  Q5(t)  (flux  pulse) 

7T2(i2»t)  *  0  back-side  heat  loss) 

Tj(0,t)  ■  T2(0,t)  (no  interfacial  resistance) 
K,7T,(0,t)  -  K,7T,(0,t) 


(13) 

(14) 


(15) 

(16) 
(17) 


(zero  divergence  at  Interface) 


(18) 


Notice  particularly  that  the  temperature  pulse  in  the  IC  of  the  Sturm- 
Liouvllle  problem  has  been  replaced  by  a  flux  pulse  in  the  BC  of  this  Laplace 
problem.  Combining  the  above  equations,  we  have  the  following  results: 

PDE  &  IC: 


Kj  — 2 - {sTj  -  0}  «  0 

dx'^ 


K2  — 2 - ^2^2  ^*"^2  “  0}  -  0 

dx 


—  T  =0 

“1  ‘ 


1  I  dx 


s  >  0 


X  =  -i. 


X  » 


T,  (0,s)  -  T„  (0,s) 


s  >  0 


X  ■  0 ,  s  >  0 


Letting  cOj.  =  /s/a^,  the  DDEs  in  Eqs.  (21)  and  (22)  have  the  solutions 
cosh  u>^x  +  slnh  u^x.  Using  the  BCs  in  Eqs.  (23)  through  (26),  one 


obtains 


where  the  matrix  a  has  the  elements  a^  (i,j  «  1,2, 3, 4),  and  3^3  *  -sinh  ^^2* 

aj^  -  -cosh  (^^12  ^21  ~  *23  “  *  *32  “  *34  ”  “^^2“2*  ®41  *  “l^l’ 

*42  “  “Cosh  all  other  a^^j  being  zero.  Furthermore,  p  *  Q/K^Uj^. 

Solving  by  determinants, 


det  (a)  *  ^2*^2  “2^2 

+K^a)j^  sinh  cosh  <^>2^2 


Aj  =  pK^u)^  cosh  (1)2 (i) 
*  -PK2OJ2  sinh  0)212/^®*^ 

^2  “  cosh  012*^2^*^®^ 

B2  *  -pK^iOj^  sinh  012^2/^®^ 


Introducing 


l^2^2^2  / - 

- -  and  /U.s  =  1.  u. 

K^PjCj  i  i  i 


and  noting  that  through  (32)  become 

det  (a)  ■  /s  /K^PjCj  [sinh  /U^s  cosh  /ujs  +  a  cosh  /U^s  sinh  /U^]  (33) 


^  *  •  «  •  -  A  ^•.  •  A  •  ■- 


(34) 


Hence 


Q  cosh  (a) 

Q  (-  o  sinh  (a) 

Q  cosh  (a) 

Q  (-  sinh  (a) 


"  det^(a)  “2^2  ~  °  “2^2 


^2  “  Te~t^(a)  “*2^2  “2^^  ~  “2^2  “2*^  (39a) 


'^2  “  det^(a)  ’  ’'^'^2^^ 


(39  b) 


Combining  Eqs.  (33)  and  (39b)  and  introducing  T^  = 


,  we  have 


^PlCi 


cosh  [/U„s  (1  -  x/&_) ] 


T  1/2  —  - 

0  s  sinh  /U^s  cosh  /U2S  +  a  cosh  /U^s  sinh  /U2S 


On  the  entire  complex  plane  the  series  expansions  of  the  factors  in  the 
denominator  vary  as  /s  (/s  (1  +  even  powers)]  (1  +  even  powers).  Therefore 
the  denominator  behaves  as  s  near  s  *  0  (simple  pole  at  s  >  0).  The  nature  of 
the  remaining  roots  of  the  denominator  has  been  investigated  (pp.  324-326  of 
Ref.  14,  pp.  409-410  of  Ref.  15  and  Appendix  A  of  Ref  16).  There  are  no 
complex  roots.  The  function  T2/TQ  is  meroraorphic,  with  poles  at  s  ■  0  and 
®  *  ®n»  ”  “  1»2»3 . 

The  Mlttag-Lef f ler  Theorem  states  that  any  meromorphic  function  can  be 
expressed  as  a  sum  of  an  entire  function  and  a  series  of  rational  functions. 

In  our  case,  letting  f(s)  ■  '^2'^^0 


r/_\  _  ?  n  ^  1  f  f(s')ds' 

+  2^1  J  S-  -  S 
n-0  n  c 


where  is  the  residue  of  the  n^^/pole.  The  contour  Integral  in  Eq.  (Ala) 
can  be  shown  to  be  zero  as  in  standard  texts  on  complex  analysis.  Thus 

^  i. 


£(.).- 11.  st(s)+J  f- -(V) 


where 


cosh  [/U^s  (1  -  x/jt^)] 


Y<s)  =  sinh  /ujT  cosh  /U^  +  a  cosh  /uJT  sinh  /ujs  (A3 


Thus,  letting  X  *  obtain 


'^2  1  •  X  (s^) 

’’^0  s  /U^  (X  +  a)  n-1  (s  -  s^) 


Using  the  inverse  Laplace  operator, ,  one  obtains 


T-  ,  »  X  (s  )  s  t 

^  (X  *  o)  *„.^1  ‘ 


Consider  the  following  definitions: 


•>  -j.'- V  -  .-j.  • .  -j.  -  !."■ 


(46  b) 


Xo  +  1 
X  +  1 


o 


H 

X 


(46  c) 


U,  -  — ,  1  -  1,2 
1 


(46d) 


After  considerable  algebra.  It  results  that 


AjPlCi  + 


1  +  2  I 


-B^t/U 

e  cos  [B  (1  -  x/Jt  )] 

n  2 


cosXS  cosB  ~  JlsinXB  slnB 
n-l  n  n  n  n 


(47a) 


V(t) 


T2  (^2**:) 


2  ,max 


(47  b) 


where  T2^^3jj  •  Q/(AjP^Cj  +  A2P2^2  equation  that  gives  the  roots  of  the 

denominator  in  Eq.  (40),  when  recast  in  terms  of  circular  functions,  deter¬ 
mines  the  real,  positive  eigenvalues  0  “  “  ’^~^2®n 

sin  XB  cosB  +  o  cos  XB  sinB  “  0  (47c) 


Furthermore,  when  the  slabs  have  identical  parameters. 


T  ■ 


1  +  2  [  (-1)" 

n»l 


-n^x^Ojt 

At? 


cos 


(1 


x/  ^2) 


(47d) 


which  corresponds  to  the  V  result  given  in  the  homogeneous  sample  analysis. 
This  shows  that  the  flux-pulse,  dual-slab,  Laplace  analysis  can  reduce  to 
the  temperature-pulse,  homogeneous-slab,  Sturra-Llouvllle  analysis  for  an 


appropriate  choice  of  parameters.  Finally,  if  V(t)  or  is  measured  and 

compared  to  Eqs.  (A7a)  or  (A7b),  one  unknown  material  parameter  can  be 
obtained  by  standard  computer  techniques. 
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I  -  ^  ^  \ 


>  Aa-Xa 


A  . 


V.- i.~-  c.  c-  .•uV  .*■ 


III.  ANALYSIS  OF  THE  RADIANT-ENERGY  PENETRATION  EFFECT 


Suppose  that  the  Incident  electromagnetic  energy  Is  converted  to  heat 
throughout  the  depth  of  the  first  layer  of  a  two-layer  composite  according  to 


-k(x  +  tj) 
c  e 

-k£  ’ 

(1  -  e  M 


-  <  X  <  0 


The  normalization  was  chosen  to  ensure  that  Q  energy  units  per  area  are 
Incident.  This  Is  the  same  Q  as  In  the  standard  solution  above,  In  which  all 
the  electromagnetic  energy  is  converted  to  heat  at  the  front  interface.  We 
will  first  analyze  what  effect  this  penetration  will  have  on  the  observable 
quantity  tjy2* 

Equations  (11)  through  (18)  are  the  same  as  before,  with  the  following 
exceptions: 


Tl(x,0) 


-k(x  +  £j) 


-k£  ’ 

PjC^d  -  e  ) 


-  £j  <  X  <  0 


-K^7T^(-£j,t)  -  0 


Now  Eqs.  (19)  through  (26)  are  the  same  as  before,  with  the  following 
exceptions : 


j2-  -k(x  +  £, ) 

K  ill  -  p  c  -  s T  -  ^ 

*^1.2  ^11  ^1  -k£,. 


The  solution  to  the  ODE,  Eq.  (19''),  is  somewhat  different  than  before: 


Tj(x,s)  -  cosh  ujj^x  +  sinh  oj^x  + 


Qke 


-k(x  +  £j) 


“k^i  2  2 

K^(l  -  e  (a)j  -  k^) 


Now  the  matrix  equation,  Eq.  (27)  becomes 

Ai\  1°  ] 


All  the  entries  of  a  are  the  same,  but  now 


-k£, 


-Qke 


Kj  (1 


-kJl 


-  e 


f  2 


k^) 


-k£ 


4Qk  e 


1 


“^■^1  2  2 
(1  -  e  (u‘  -  k^) 


-Qk 


2  2 

(1  -  e  (o)^  -  k^) 


-k£. 


We  will  assume  only  a  limited  penetration  consistent  with  e  <  0.05, 


(19'-) 

(27') 

(A9a) 

(49b) 

(49  c) 

or  k£,  >  3 


To  a  reasonable  approximation,  therefore,  we  set  e  ,  p",  and  p"  equal 
to  zero ,  and 


In  the  limit  that  k  +  »  (delta  function,  or  surface,  energy  deposition), 

p  +  Q/KjU)j  as  In  the  standard  solution  above.  Therefore,  we  can  follow  the 

changing  solution  as  surface  deposition  Is  continuously  changed  to  deposition 

in  depth.  We  can  do  this  rigorously  if  nonzero  p"  and  p'  are  considered. 

Here  we  consider  only  the  case  of  small  penetration,  kjt,  >  3. 

1  ~ 

Since  the  problem  has  only  changed  to  the  extent  of  a  modification  of  the 
factor  Q  in  the  quantity  p.  It  is  seen  Immediately  that 

T2  cosh  [/U-s  (1  -  x/i  )] 


0  s  '  (1  -  s/k  Oj)  [sinh  /U^s  cosh  /u^  +  o  cosh  /UjS  sinh  /U^s] 

Hence  the  change  In  the  problem  which  Is  due  to  finite  radiant-energy  penetra 
tlon  Is  seen  In  a  change  in  the  pole  structure,  as  shown  in  Fig.  1. 

Thus  the  previous  solution,  Eq.  (Alb),  Is  modified  to  become 


f(s)  1 


_o+  y  — + 

s  s  -  s  s  -  s, 

n«l  n  k 


Since  k)l^  >  3,  k  is  positive,  and  Sj^  >  0.  Now,  considering  Eqs.  (A2),  (43), 
and  (40"), 


°  /V2  (X  +  a) 


(51b) 


r 

n 


-^k  » 
<*n  - 


X  <\) 
'k  ■  ■  ®k  Y  (s^) 


(51c) 


-kt, 

Wej^gote  that  tj^  is  proportional  to  e  ♦0,  since  we  have  taken 
e  ^  =  0  as  an  approximation.  Thus  the  j^^st  term  in  Eq.  (Alb")  must  be 
ignored.  Only  if  p"*0,  p^^O,  and  e  ^  ^  0  in  Eq.  (27 ')  would  it  be  con¬ 
sistent  to  keep  such  terms.  In  our  approximation,  then,  Eq.  (41b'’)  differs 
from  Eq.  (41b)  only  through  the  ~Sj^/(Sjj  -  Sj^)  multiplicative  factor  in 
Eq.  (51b). 

=  1,2 


Using  Sj^  =  k^a^  and  s^^  »  becomes 


2 ,  max 


1  +  2  I 


1 


n-1  1  + 

n  2  1 


2 

cos  [6^  (1  -  x/jl^)] 

cosXB  cosB  -  n  sinXB  sinB 
n  n  n  n 


(52) 


Equation  (52)  is  the  master  equation  for  the  p"  =  p'  =  0  approximation  to  the 
finite-energy  penetration  problem.  Two  tasks  now  remain:  First  it  will  be 
demonstrated  that  this  energy  penetration  can  significantly  affect  experi¬ 
mental  results.  Then  it  will  be  indicated  how  the  computer  data-reduction 
algorithm  might  be  modified  to  include  codetermination  of  k  along  with  the 
usual  parameter  evaluation. 

The  significance  of  Eq.  (52)  can  be  judged  by  evaluating  it  for  the  case 
of  identical-layer  parameters  and  in  the  first-term  (n  =»  1)  approximation 
analogous  to  Eq.  (9) : 


1  .  2  e~'"^°l^l/2'^^^l 

2  “  .  .  2,,, 2.2 


2 

n  a 


1  +  ir‘/4k''Jl 


1 


4£,  2 

—  [in4  -  £n(l  + 


2  2 

4kjr 


)] 


’1/2 


23 


(53a) 


(53b) 


Letting  rise  time  for  the  spatial  delta-function  energy  deposition 

[  Eq .  ( 9 )  ]  ,  we  have 


Tl 

2  2  ^ 


(53c) 


Letting  SF  be  the  fractional  deviation  (shortfall)  of  the  rise  time  in  the 
finite-deposition  case,  we  obtain 


t  -  t^0> 
1/2  *^1/2  _ 

t(0) 

1/2 


And  +  Ti^/4k^Aj) 


(53c) 


for  the  approximation  kA^  >  3.  Figure  2  shows  the  result  for  the  shortfall 
of  the  rise  time  versus  kA^  or  the  characteristic  penetration  depth  1/k. 

The  approximations  used  in  the_g2*'®t*’'^ction  of  Fig.  2  could  be  removed 
by  the  rigorous  (p'"  ^0,p'*0,e  ^#0)  solution  of  Eq.  ( 27 ") ,  the 

subsequent  Inclusion  of  the  rj^  residue  in  Eq.  (41b'’),  and  the  use  of  the 
entire  series  in  Eq.  (52).  However,  even  with  those  refinements  the 
-SF  (kA^  =  3)  is  still  expected  to  be  in  the  10  to  20%  range,  even  though  the 
energy  reaching  the  midplane  of  the  sample  is  only  e”^  =  5%  of  that  converted 
to  heat  at  the  front  interface. 

We  see  that  the  finite  penetration  of  radiation  may  cause  a  significant 
rise-time  shortfall.  How  can  this  effect  be  allowed  for  in  data  reduction? 

In  one  standard  form  (p.  493  of  Ref.  10)  of  the  flash-technique  data 
reduction,  the  eigenvalue  equation  [Eq.  (47c)]  is  solved  with  the  measured 
auxiliary  parameters  and  a  first  guess  at  a, ,  which  is  chosen  as  the  unknown. 
The  first  guess  is  “  0.1388  ^i/^j/2’  *  (Presumably  •  0.1404 

fj/ti/2  [Eq.  (9)1  would  also  lead  to  a  convergent  series.)  The  deviation  from 
V(t  *  ^1/2^  calculated,  and  an  Iteration  algorithm  is  devised  which 
terminates  when  a|  ^  ^  ■  y  ( oj  +  a|  ^)  converges  to  a  reasonable  limit. 
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The  Magnitude  of  the  Percent  Shortfall  (-SF)  of  the 
Measured  Rise  Time  in  the  Thermal  Diffusivity  Measure¬ 
ment  of  a  Homogeneous  Layer  of  Width  2i. ,  Graphed  as  a 
Function  of  ki  and  of  the  Characteristic  Penetration 

^  ....  -ki 

Depth  of  1/k.  In  our  approximation  (p^'-p'^e  !■ 

kJl,  >  3,  and  only  the  n  ■  1  rise  times  are  compared. 


Bulmer  and  Taylor  use  only  the  ^1/2  experimental  output,  because  they  feel 
that  curve  fitting  to  more  points  requires  using  extremes  of  the  experimental 
curve  which  are  inherently  less  accurate. 

We  would  argue  that  the  entire  experimental  V(t)  curve  contains  some 
relevant  information,  and  that  at  the  very  least  one  additional  parameter,  k, 
should  be  derivable  from  this  information.  We  propose  that  a  figure  of  merit 
RMS^  (k-»-<»,  ^),  be  formed  for  the  first  determination,  aj  ^  ^ ,  in  the 

following  manner: 


(k  .  aj  *  ')  -i  I  |V  (t__)  -  V^(t__)r 


where  N  points  t^^  on  the  time  scale  are  chosen.  Vgj^p(t)  is  the  experimental 
trace,  and  Vj^(t)  is  given  in  Eq.  (52),  in  this  case  with  k  ■♦•  ».  Now  a 
different  initial  guess  can  be  made  for  using  a  finite  value  of  k: 

41^  2 

a,  ■  —  [kn4  -  4n(l  +  — ^)  ]  (55 

’  '^1/2  *1 

1“*  +  1  ?  f*  +  1 

After  convergence  to  a  new  a.  ,  the  quantity  RMS  (k,  a,  )  can  be 

2  1'  +  1  ^ 

calculated.  When  RMS  (k,  )  reaches  a  minimum  as  a  function  of  k,  that 

is  the  best  k  for  the  measurement. 


IV.  SUMMARY  AND  CONCLUSIONS 


Current  and  future  aerospace  and  construction  Industry  materials  require 
efficient  thermal  analysis  techniques.  Including  the  flash  technique  for 
thermal  dlffuslvlty.  This  Is  not  possible  unless  the  data  can  be  adequately 
reduced  for  polymer  composites  and  porous  blankets.  This  paper  has  been 
concerned  with  radiant-energy  penetration  Into  a  layered  sample  subjected  to 
flash-technique  conditions.  A  mathematical  modification  of  the  standard 
flash-technique  solution  for  layered  structures  was  presented  which  allows  for 
an  exponentially  decaying  penetration  of  radiant  energy  Into  the  first  of  two 
layers.  The  following  conclusions  were  drawn: 


(1) 


(2) 


A  new  (simple)  pole  In  the  Laplace-transformed  solution  appears 
on  the  complex  s  or  S  plane.  Thus  the  Inverse  Laplace  transform, 
which  gives  the  response  function  (back-face  temperature  versus 
time),  contains  an  ^ddlt^onal^term.  Furthermore,  there  Is  a  small 
modification  (1  +  6^/U2k"‘oj)  to  the  previous  terms. 

The  appearance  of  a  new  pole  signifies  In  general  the  emergence  of 
a  new  kind  of  physical  behavior,  because  the  residue  of  the  new 
pole  may  in  general  be  significantly  larger  than  those  of  many  of 
the  original  poles.  However,  to  facilitate  the  evaluation  of  the 
Importance  of  thjLj  effect,  a  small-penetration  approximation  was 
made  In  which  e  was  set  equal  to  zero. 


This  Is  approximately  true  for  k£j  >  3  (e  1  <  0.05)  .  In  this 
approximation,  t^e  new  term  Is  suppressed  and  only  the 
(1  +  S  /U-k  a.)  modification  retained.  However,  It  must  be  kept 
in  mind  that  tor  more  extensive  radiant  energy  penetration,  this 
approximation  would  not  apply,  in  which  case  the  new  term  In  the 
solution  could  specify  a  radical  departure  of  the  response 
function  from  Its  previous  behavior.  Furthermore,  extensive 
penetration  might  also  require  the  exponential  function  Q(x)  to 
extend  over  more  than  one  layer. 


(3)  A  numerical  example  was  evaluated  for  the  e  1  ~  0  approximation 
in  the  case  of  Identical  layer  parameters  and  the  first-term 
(n  ■  1)  approximation  to  the  back-face  rise  time.  It  was  found 
that  the  rise  time  shortfall  (l/k  ■  33%  of  %j)  •>  18%  of  the  rise 
time  that  would  obtain  if  all  the  radiant  energy  were  deposited  on 
the  surface.  Here,  l/k  Is  the  characteristic  penetration  depth  of 
the  exponential  radiant-energy  deposition  function.  Thus,  even 
when  the  radiant  energy  deposited  at  the  sample  midplane  Is  only 
5%  of  that  deposited  on  the  surface,  an  18%  correction  is 
necessary. 


The  standard  fJa^hYtechnique  data-reduction  algorithm  considers 
Iterations,  aj  - -j  [oij  +  '],  which  are  required  to 

converge  to  a^ reasonable  limit.  Om.y  the  tj^^2  experimental  data 
are  used.  If  the  entire  V(t)  response  function  were  used,  how¬ 
ever,  a  figure  of  merit  might  be  defined  that  allows  for  the 
experimental  determination  of  k. 

Modifications  such  as  the  ones  reported  here  are  expected  to  open 
up  the  flash  technique  to  the  analysis  of  advanced  polymer 
composite  and  insulating  blanket  material. 
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LIST  OF  SYMBOLS 


Aj  , 


a 


c 

f(s) 

H 


constant  coefficients  In  t(ie  solution  of  the  ordinary 
differential  equation  for  T2. 

matrix  of  the  coefficients  B^  arising  from  the  boundary 
conditions  (see  Eq.  (27)] 

specific  heat 

-  T^Cx.s)/!^ 

=«  X/o 


i 


layer  index  for  the  composite  slab  (layer  1  is  hit  by  the 
incident  flash) 


J 


thermal  flux 


K 

1/k 


i£  , 


L 

"1 

Pt  P  I  P 


Q 


^n 

RMS^(k,  aj  *  S 


s 


s 


n 


thermal  conductivity 

characteristic  penetration  depth  of  the  radiant  energy  into 
the  first  layer  of  a  composite  slab  (k  +  «*  for  surface 
deposition  only) 

Laplace  (inverse  Laplace)  operator 
thickness  of  the  (homogeneous)  layer 
thickness  of  the  i^^  layer 
nonhomogeneous  entries  in  Eqs.  (17)  or  (17“’) 
heat  per  unit  area  deposited  on  sample 
residue  of  f(s)  at  s  *  0 

residue  of  f(s)  at  s  ■  Sjj  (n  »  1 ,2 ,3 ,  ...) 
residue  of  f(s)  at  s  Sj^ 

flgure^o^  ^erlt  for  deviation  of  Vy^(t)  from  Vgjjp(t),  using 

“1  “  “1 

Laplace-conjugate  variable  to  the  time  t 

location  of  the  simple  poles  of  T2/TQ  on  the  complex  s-plane 
(n  “  1,2,3,  ...)  (Sq  «  0  simple  pole  evaluated  separately) 
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LIST  OF  SYMBOLS  (Continued) 


SF 


Ti(x,t) 


Tj^(s,t) 


■1/2 


^2 ,max 


.(0) 


1/2 


X(8),Y(s) 


*  ■  the  new  simple  pole  that  appears  as  a  result  of 

finite  radiant-energy  penetration 


fractional  shortfall  of  the  measured  rise  time  (penetration 
case)  below  that  occurring  for  surface-energy  deposition  only 
[see  Eq.  (53c)] 


temperature  In  layer  1  at  position  x  and  time  t 
Laplace  transform  of 


characteristic  rise  time  of  the  back-face  temperature  of  the 
sample  In  the  flash  technique 


Q//KjPjC^ 


back-face  rise  time  when  there  is  surface-energy  deposition 
only 


- 

/u][s 

V(t) 

back-face  response 

function 

'ezpf'' 

measured  back-face 

response 

function 

back-face  response 
finite  k 

function 

[see  Eq.  (52)]  for  the  case  of 

✓Uj/U2 


useful  functions  for  residue  evaluation  [see  Eqs.  (42) 
and  (43)] 


thermal  dlffusivlty 


the  value  of  after  1  Iterations  of  the  standard  data- 
reductlon  algorithm  for  the  flash  technique 


eigenvalues  of  the  boundary-condition  problem  (■  /-uTs  , 
n  -  1,  2,  3,  ...)  " 


differential  operator 
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LIST  OF  SYMBOLS  (Continued) 


5 

P 

a 

0 


Dirac  delta  function 
density 
*  K2U)2/KjU1  - 
(Xo  +  1)/(X  +  1) 

=  /s/a^ 


Vl'VO.-Vi.'.- 


LABORATORY  OfERATIOWS 


The  lAboretory  Operetlone  of  The  Aeroepace  Corporation  la  conducting 
cxperlaental  and  theoretical  Inveatlgatlona  neeeaaary  for  the  evaluation  and 
application  of  scientific  advances  to  neu  ellltary  space  systees.  Versatility 
and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  person¬ 
nel  In  dealing  with  the  nany  problees  encountered  In  the  nation's  rapidly 
developing  space  systees.  Expertise  in  the  latest  scientific  developecnts  Is 
vital  to  the  accoapllahaent  of  tasks  related  to  these  problees.  Tha  labora¬ 
tories  chat  contribute  to  this  research  arc: 


Aerophysles  Laboratory:  Launch  vehicle  and  reentry  fluid  aechanlcs,  heat 
transfer  and  flight  dynaalca;  chealcal  and  electric  propulsion,  propellant 
chealatry,  envlronaental  haxards,  trace  detection;  spacecraft  structural 
aechanlcs,  contaalnatton,  thsraal  and  structural  control;  high  caaperature 
theraoaechanlcs,  gas  kinetics  and  radiation;  cw  and  pulsed  laser  dcvelopaent 
Including  chealcal  kinetics,  apeccroacopy.  optical  resonators,  bean  control, 
ataospherlc  propagation,  later  effects  and  counteraeaaures. 


Chealatry  and  Physics  Laboratory:  Ataospherlc  chealcal  reactions,  atao- 
tpherlc  optics,  light  scattering,  state-specific  chealcal  reactions  and  radia¬ 
tion  transport  In  rocket  pluses,  applied  laaer  spectroscopy,  laser  chemistry, 
later  optoelectronics,  solar  cell  phyalct,  battery  elcctrochealttry,  space 
vacuum  and  radiation  effects  on  aatertals,  lubrication  and  turfsee  phenoaena, 
themlonlc  emission,  photosensitive  atterlals  and  detectors,  atoalc  frequency 
standards,  and  envlromaental  chealatry. 


Electronics  Research  Laboratory:  Wcroclectronies,  GaAs  low  noise  and 
power  devices,  sealconductor  lasers,  electroaagnetlc  and  optical  propagation 
phenoaena,  quantum  electronics,  laser  coaaunl cat Ions,  Ildar,  and  electro- 
optics;  coamunlcatlon  sciences,  applied  electronics,  sealconductor  crystal  and 
device  physics,  radiometric  laaglng;  mllllaeter  wave,  microwave  technology, 
and  RP  systeas  research. 


Information  Sciences  Research  Office;  Prograa  verification,  prograa 
translation,  perforaance-sensltlve  system  design,  distributed  architectures 
for  spaceborne  computers,  fault-tolerant  computer  systeas,  artificial  Intel¬ 
ligence  and  aicroclectronlcs  applications. 


Haterlals  Sciences  Laboratory:  Dcvelopaent  of  new  aaterlals:  aetal 
aatrlx  coaposltes,  polyaers,  and  new  foras  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  tellablllty;  fracture  aechanlcs  and 
stress  corrosion;  analysis  and  evaluation  of  aaterlals  at  cryogenic  and 
elevated  ceaperatures  as  well  as  in  space  and  eneay-lnduced  environments. 


Space  Sciences  Laboratory:  Msgnetospherlc,  auroral  and  cosalc  ray  phys¬ 
ics,  wsve-partlcle  Intersctlons,  asgnetospherlc  plasas  waves;  ataospherlc  and 
Ionospheric  physics,  density  and  coaposltlon  of  the  upper  ataosphere,  remote 
sensli^  using  atmospheric  radiation;  solar  physics.  Infrared  astronoay. 
Infrared  signature  analysis;  effects  of  solsr  activity,  aagnetlc  atoms  and 
nuclear  explosions  on  the  earth's  ataosphere.  Ionosphere  end  aagnetosphere; 
effects  of  elaccromagnetlc  and  particulate  radiations  on  space  systeas;  space 
Inatruaentatlon. 
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